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The transition from a normal to cancerous cell requires a number of highly specific mutations that 
affect cell cycle regulation, apoptosis, differentiation, and many other cell functions. One hallmark of 
cancerous genomes is genomic instability, with mutation rates far greater than those of normal cells. 
In microsatellite instability (MIN tumors), these are often caused by damage to mismatch repair 
genes, allowing further mutation of the genome and tumor progression. These mutation rates may 
lie near the error cateistrophe found in the quasispecies model of adaptive RNA genomes, suggesting 
that further increasing mutation rates will destroy cancerous genomes. However, recent results have 
demonstrated that DNA genomes exhibit an error threshold at mutation rates far lower than their 
conservative counterparts. Furthermore, while the maximum viable mutation rate in conservative 
systems increases indefinitely with increasing master sequence fitness, the semiconservative threshold 
plateaus at a relatively low value. This implies a paradox, wherein inaccessible mutation rates are 
found in viable tumor cells. In this paper, we address this paradox, demonstrating an isomorphism 
between the conservatively replicating (RNA) quasispecies model and the semiconservative (DNA) 
model with post-methylation DNA repair mechanisms impaired. Thus, as DNA repair becomes 
inactivated, the maximum viable mutation rate increases smoothly to that of a conservatively repli- 
cating system on a transformed landscape, with an upper bound that is dependent on replication 
rates. On a specific single fitness peak landscape, the repair-free semiconservative system is shown 
to mimic a conservative system exactly . We postulate that inactivation of post-methylation repair 
mechanisms are fundamental to the progression of a tumor cell and hence these mechanisms act as 
a method for prevention and destruction of cancerous genomes. 

PACS numbers: 87.14Gg, 87.17.-d, 87.23.-n 



I. INTRODUCTION 

Cancer has presented itself as one of the most difficult 
challenges science has ever faced. The complexity of the 
disease, experimental obstacles, and the vast array of tu- 
mor types have made characterization of the many facets 
of tumor progression a slow process. It is now understood 
that this progression requires the alteration of numerous 
genes, as a genome progresses from its normal state to a 
full-blown cancer celli. One important aspect of the can- 
cerous genome lies in its genetic instability. All cancerous 
genomes display either high mutation rates (in MIN tu- 
mors) or chromosomal instability (in CIN tumors 

One of the most successful theoretical methods for 
studying genomic evolution at high mutation rates has 
been Eigen's quasispecies model™. This model considers 
an explicit population of genomes, each made up of L 
nucleotides chosen from an alphabet of size S, usually 
chosen to be two for simplicity or four to model the nu- 
cleotides in nature. These genomes replicate, mutate and 
compete on a chosen fitness landscape, a unique mapping 
of genotype to fitness. This is often accomplished by as- 
signing different replication rates to each possible genome 
and setting all death rates to be equal. The model has 
yielded a number of impressive and experimentally ver- 
ified predictions^iSiSii and has recently been used as the 
basis for novel anti-viral therapies^. The main predic- 
tion lies in the idea of an error catastrophe. Below a 
threshold mutation rate, dubbed the "error threshold" , 
the population evolves, independent of starting condi- 
tions, to a distribution of genomes near the sequence of 



maximal fitness, often called the master sequence. Above 
the threshold mutation rate or the "error threshold" , the 
population reaches a random distribution with no dis- 
cernible master. This crossover is depicted in Fig. 1. 
These ideas provide a method for destroying RNA-based 
viral genomes. Viruses are expected to evolve a mutation 
rate slightly below the error threshold so as to main- 
tain the capacity to rapidly adapt without surpassing 
the error threshold and becoming inviable. Hence, by in- 
creasing the mutation rate of the species, the virus can 
be destroyed, and this technique has been successfully 
applied^. 

These ideas have recently been suggested to apply to 
cancer cells-. Cancer and RNA viruses share genetic 
instability in the sense that both are rapidly mutating 
and recent work has focused on the idea that mutagens 
may push cancer cells past the error threshold in a sim- 
ilar manner. Support for the idea that the quasispecies 
model can be applied to complex cellular genomes comes 
from recent studies that yielded accurate qualitative and 
quantitative predictions for complex systems such as the 
adaptive immune system.^. 

However, past work on the quasispecies model has 
focused on conservatively replicating systems such as 
RNA-^^i-^-^i-^^i-^'^i-^'^i-^^. In these systems, single stranded 
genomes are copied to produce a new, possibly error 
prone, strand without affecting the original. In semicon- 
servative systems like DNA, double stranded genomes un- 
zip to produce two single strands, each of which is copied 
to produce a new complementary strand by Watson- 
Crick base-pairing. A variety of mismatch repair enzymes 
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FIG. 1; A schematic diagram illustrating the error catastro- 
phe predicted by the quasispecies model. The concentration 
of sequences of Hamming distance i from the master sequence 
is represented by Xi and plotted against the mutation rate, e. 
At low concentrations, the master sequence (of Hamming dis- 
tance 0), dominates the population, but is surrounded by a 
cloud of closely related genomes . Above the error threshold, 
this clustering disappears, and we see a random distribution 
of genomes, where each Hamming class has a concentration 
proportional to its size. 
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then repair any errors in the new strand, keeping the ef- 
fective error rate low. A few errors remain, though, and 
these, as well as extrinsic mutations induced by UV radia- 
tion or other mutagens, are repaired by post-methylation 
repair enzymes that cannot distinguish between the new 
and old strands. Thus, some of these base pair mis- 
matches are repaired in the old strand and the original 
strand is not conserved-- as shown schematically in Fig. 
2. 

The quasispecies model has recently been extended to 
incorporate this behavior—. It was found that, contrary 
to popular assumption, the semiconservative system dis- 
plays fundamentally different behavior than the conser- 
vative system. In particular, on a single fitness peak land- 
scape, the semiconservative error catastrophe occurs at 
significantly lower mutation rates than the conservative 
case-"^^. Particularly interesting is that, in the conserva- 
tive case, the maximum viable mutation rate increases 
without bounds with increasing master sequence repli- 
cation rate, while the semiconservative system reaches a 
threshold valueii, and this has been confirmed by sim- 
ulation for finite populations sizes and genome lengths. 
Thus, for the conservative case, it is always possible to 
"out-replicate" the error threshold. That is, for any given 
mutation rate, there exists a relative fitness for the mas- 
ter sequence such that, if the master sequence has that or 
greater fitness, the error catastrophe is avoided. This is 
not true for semiconservative genomes, where there ex- 
ist mutation rates that cause the error catastrophe for 



FIG. 2: A schematic model of (a) conservative and (b) semi- 
conservative replication. Non-methylated strands are bolded 
and errors are circled. 



any value of the master sequence fitness. We note in 
passing that a conservative system will also follow this 
behavior if the possibility of unrepaired extrinsic muta- 
tion is incorporated. The fundamentally different nature 
of the semiconservative error catastrophe has numerous 
implications^, but is particularly pertinent to the study 
of cancer. 

Within the conservative paradigm, it is reasonable to 
assume that cancer cells are capable of maintaining a vi- 
able population of rapidly mutating genomes, as rapid 
replication rates are one of the hallmarks of cancerous 
cellai. However, the recent results on the semiconserva- 
tive system present a paradox. The mutation rates in 
MIN cancer cells, known to be 50-1000 times higher than 
those of normal cellsiSiiS, certainly lie higher than any 
reasonable value for the low semiconservative threshold 
(for example, the single fitness peak landscape yields 1.39 
errors/genome/replication as an upper bound for the er- 
ror threshold in a long semiconservative genome, while 
cancer cells display error rates over three orders of mag- 
nitude greater). Furthermore, the rapid replication rates 
that allow such high mutation rates in the conservative 
case provide no help, as the maximum allowed muta- 
tion rate cannot exceed a rather low threshold value, no 
matter how fast the cancer cells replicate. Hence, these 
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recent results appear to present a paradox: rapidly mu- 
tating genomes are prevalent in cancerous cells, but such 
high mutation rates should exceed the error threshold 
and hence yield inviable genomes. 

In this paper, we address this paradox and demon- 
strate that a scmiconscrvativc system can mimic a con- 
servative population through the degradation of post- 
methylation lesion repair. In section II we discuss the 
isomorphism between conservative replication and semi- 
conservative replication without lesion repair. In section 

III we look at the implications of this result and in section 

IV we present our conclusions. 



II. DNA REPAIR AND SEMICONSERVATIVE 
REPLICATION 

As discussed above, DNA replication can be considered 
a three part process; unzipping, complementary strand 
creation and mismatch repair. Afterwards, any remain- 
ing mismatches, as well as damage caused by environ- 
mental conditions, are repaired by a set of repair en- 
zymes. Global genomic repair (GGR) fixes lesions, errors 
and mismatches along the entire genome, while transcrip- 
tion coupled repair (TOR) subjects the expressed portion 
of the genome to more careful scrutiny and repair. 

In appendix A, we use the quasispecies equations 
to demonstrate a mathematical isomorphism between 
a population of conservatively replicating genomes and 
semiconservatively replicating genomes without any le- 
sion repair. In this case, the semiconservative system 
behaves, in essence, like a conservatively replicating sys- 
tem on a transformed landscape. Each single stranded 
genome produces one, possibly error prone, complemen- 
tary copy. Although mismatch repair may keep the effec- 
tive error rate low, the lack of lesion repair ensures that 
the original strand is unaffected by these errors. Hence, 
each genome replicates, in essence, conservatively, but 
with the added wrinkle that each single stranded genome 
remains attached to the strand that either created it or 
that it most recently created, yielding a system that repli- 
cates conservatively on a transformed fitness landscape. 

To make this more rigorous, appendix B presents the 
full solution to the semiconservative quasispecies evolv- 
ing on a specific single fitness peak landscape. This is 
plotted in Fig. 3, along with the conservative and semi- 
conservative solutions to the same problem. While the 
semiconservative error threshold clearly plateaus at high 
cr, the repair-free semiconservative case mimics a conser- 
vative system, as the error threshold increases indefinitely 
with increasing a. 

Lastly, it is important to consider the case where le- 
sion repair is partially active, as complete degradation 
of lesion repair is not likely to occur in nature. In ap- 
pendix C, the single fitness peak quasispecies is recon- 
sidered, this time with partially active lesion repair. The 
error threshold is shown to increase smoothly from the 
semiconservative to the conservative threshold as shown 
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FIG. 3: The value of the error threshold (eo) vs. the fitness of 
the master sequence relative to the rest of the population (cr) 
on a single fitness peak landscape. The genome length is set 
to AT = 1 X 10*. Conservative, semiconservative and semicon- 
servative systems without post-methylation lesion repair are 
all shown. 

in Fig. 4. This turns out to be important, and will be 
discussed in the next section. 
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FIG. 4: The error threshold (eo) for a semiconservatively 
replicating population on a single fitness peak landscape vs. 
the probability of post-methylation DNA repair (A). Here, 
= 1 X 10* and (7 = 1 X 10^ 



III. DISCUSSION 

As discussed above, recent results on the semiconser- 
vative quasispecies has presented a paradox in tumor 
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progression. In the last section, we presented a possi- 
ble resolution. Although extremely high mutation rates 
and genetic instability are found in all cancer cells, these 
mutation rates cannot be handled by a semiconservative 
genome. However, as lesion repair begins to fail, the er- 
ror threshold increases. These increased mutation rates 
may in turn lead to further failure of the lesion repair sys- 
tem and a higher threshold, creating a positive feedback 
cycle. 

The concept of a mutator phenotype in cancer has a 
long and established historjiSS. Failure to prevent and 
repair mutations is well documented in cancer cells2i, 
yet it remains unclear to what extent each repair en- 
zyme is or is not active in any given cell. Thus, it is 
difficult to say with conviction that lesion repair fail- 
ure is indeed a prerequisite for the sustenance of ge- 
netic instability in MIN genomes in nature. However, 
there are a number of encouraging signs that this is, in 
fact, the case. Many human tumors have been found 
to be deficient in checkpoint pathways, including those 
that involve p53, pl6 and pl9ARF.^fli22,. These check- 
points are, among other things, designed to increase the 
efficiency of DNA post-methylation damage repair, and 
mice that lack these checkpoint genes display abnormally 
high levels of spontaneous tumor incidenceSS. Loss of the 
p53 tumor suppression gene has been shown to lead to 
less efficient GGR^^'^"*-^^ and mutations of the BRCAl 
gene, which enhances the GGR process, and greatly in- 
creases the risk of breast cancer in women2&. Well docu- 
mented diseases, such as Xeroderma Pigmentosumi2i42LS& 
are caused by defects in GGR and manifest themselves 
as an overwhelmingly high probability of tumor develop- 
ment. Expression profiles of pancreatic cancer cells have 
demonstrated down-regulation of DNA repair genesSS. 
All-irans-retinoic acid has been shown to prevent certain 
carcinogenic transformations by enhancing DNA repair 
through checkpoint effecta^S. Lastly, numerous human 
studies have shown positive correlations between indi- 
vidual fluctuations in DNA repair capability and cancer 
riski^. 

Although it is clear that DNA repair is linked to tumor 
suppression, we are suggesting a fundamentally different 
outlook on the problem. Rather than simply protecting 
the genome from mutations, DNA repair also prevents 
the proliferation of genomes with high mutation rates. As 
methyl-directed mismatch repair begins to fail, the error 
threshold will soon be crossed, unless DNA repair be- 
gins to fail as well. The available experimental evidence 
shows a definite correlation between repair failure and 
cancer risk, but causation is not evident. As well, alter- 
native hypotheses can explain this correlation (as repair 
failure makes genetic instability more likely, and genetic 
instability makes repair failure more probable), so more 
evidence is required. One fundamental and novel aspect 
of our hypothesis, although, of course enormously diffi- 
cult to practically implement, is that reinstating lesion 
repair in a full-blown cancer cell should lower the error 
threshold and provide the same effect as pushing the cells 



past the error threshold, without the added side effects 
associated with introducing mutagens into the body. 

One must take great care in making grand state- 
ments regarding complex biological systems from sim- 
plified physics model. Complex processes involving nu- 
merous enzymes and prefactors are incorporated into 
first-order rate constants, and the various types of DNA 
damage are ignored in favor of simple point mutations. 
Further, DNA repair covers a complex set of phenom- 
ena, rather than the simple post-methylation mismatch 
repair treated in the model. However, the quasis- 
pecies model has been impressively successful in dealing 
with a wide variety of complex systems, including finite 
populations^SiSiM, time dependent landscape g-^^i^^'^^ , 
punctuated equilibriu m^^ i^^ and even the accurate pre- 
diction of human B-cell mutation rates^ and viral 
propertiesSii. Despite its simplicity, the model seems to 
capture the robust properties of genomic evolution. Fur- 
thermore, it is successful at all mutation rates, whereas 
many theories of population genetics only work at low 
mutation rates, which obviously does not apply to genet- 
ically unstable tumor progression. Regarding the fitness 
landscape, although cancerous genomes can be highly 
heterogeneous, the single fitness peak landscape likely 
captures the general features of local behavior even on 
more complex landscapes, and can be shown to yield the 
same behavior as more delocalized landscapes. As well, 
the mathematical isomorphisms shown in the appendices 
hold for all landscapes and W matrices. 

Lastly, although the model is restricted to errors in 
the form of point mutations, these are the major source 
of genetic instability in the MIN (microsatellite insta- 
bility) tumors, which can be found in 13% of sporadic 
colon cancersSZ. and all hereditary nonpolyposis colorec- 
tal cancer. MIN tumors can display a point mutation 
rate one thousand times greater than that of a normal 
celU^i^. Other tumors display genetic instability in the 
form of GIN (chromosomal instability) with a wide vari- 
ation in chromosome number and other chromosomal 
instabilitjii2iS&. The possibility that such instability can 
be treated, to a first approximation, through the inclu- 
sion of recombination^2i^ and simulation techniques is 
the subject of future research^. Previous results in these 
areas provide reason to believe that the underlying dy- 
namics for models of GIN tumors should provide similar 
results to those obtained here. 



IV. CONCLUSIONS 

In this paper, we presented a paradox arising from re- 
cent results regarding the quasispecies model of semicon- 
servative replication. The relatively low values for the er- 
ror threshold, together with the fact that this threshold 
does not increase with increasing master sequence fitness, 
suggests that a true semiconservative system should not 
be capable of handling the exceptionally high mutation 
rates associated with cancer cells. We demonstrated that, 
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through the degradation of lesion repair, the semicon- 
servative system begins to mimic its conservative coun- 
terpart, with an increasing error threshold whose upper 
bound becomes increasingly dependent on the replication 
rate. Thus, we postulate that the failure of mismatch re- 
pair systems and the corresponding increase in mutation 
rates that are found in MIN tumors must be accompanied 
by failure of post-methylation lesion repair. Although 
we present some experimental evidence to support this, 
the simplicity of the model together with the complex- 
ity of the problem require further experimental evidence 
to fully justify our claim. Thus, we have suggested a 
new outlook to guide further experimentation and more 
complex model calculations. 

Acknowledgments 

The authors are indebted to Franziska Michor and 
Emmanuel Tannenbaum for invaluable discussions and 
a careful reading of the manuscript. This research was 
supported by an NIH postdoctoral fellowship. 

APPENDIX A: SEMICONSERVATIVE 
QUASISPECIES MODEL WITHOUT 
POST-METHYLATION DNA REPAIR 



It is important to note that, in the absence of lesion re- 
pair, (j)' is not defined by complementary base pairing 
to (j), since there is no requirement that base pair mis- 
matches be altered. We use </>') as before to de- 
scribe the probability that replication of the unzipped 
single stranded genome 4>' will produce new strand 4>. 

Hence, the quasispecies equations for a semiconserva- 
tive genome without post-methylation DNA repair can 
be written as 




Cf)a) + W{(j), (/)6))^(0a, 0fc)2:{0„,0,} - 



where we define the population fraction = 
J2(f, + ^{(t>a.<P}- l^ote that we count the 5' 3' 

and 3' 5' strand separately to avoid double counting. 
Rearranging, 

= X (/)a)^(0a,0b)a;{0„,0,} 

<P<l;<Pb 

+ J2 ^(^^ 06)^(0a>&)X{0„,0,} - f{t)x^ (A4) 



In this appendix we examine a semiconservatively 
replicating quasispecies model in the absence of any le- 
sion repair. The standard conservative model describes 
the evolution of a set of organisms, each with a genome 
4> = S1S2 ■ ■ ■ sn , where each Si represents a "letter" cho- 
sen from an alphabet of size S. The population fractions 
obey the set of differential equations'' 



~dt 



Y,A{(l)')W{cj,,cp')x^, - f{t)x^, 



(Al) 



where x^ denotes the fraction of the population with 
genome (j), A{4)) represents the fitness, or growth rate, of 
sequence (j), and W^(0, is the likelihood of creating se- 
quence 4) from (j)' by mutations. f{t) = A{if))x^ is the 
average fitness of the population, which holds the popu- 
lation size constant and introduces competition. If only 
point mutations are allowed and a genome-independent 
mutation probability e is assumed, then (j)') can be 

written in terms of the number of bases at which (j) and 
(f)' differ, the Hamming distance HD{(f>,(f)'), as 

W{c^,cj)') = (^_)^^m0')(i_g)JV-™.0'), (A2) 
»S 1 

where N represents the length of the genome. The iso- 
morphism we are about to describe holds for all W, but 
we shall limit ourselves to this manifestation in appendix 
B and C. 

For a semiconservative system, organisms are de- 
scribed by a double stranded genome {0, (f)'}, the popula- 
tion fractions as Xj^^^/j, and the growth rates as A{(j), (j)'). 



X W{(l), (j)a){A{(l)a,(j)b)X{4>^^^^} + A((/if,, 0a)a;{0(,,0„}) 



(A5) 



where the last expression is obtained by switching the 
dummy variables (pa and (pb in the second summation of 
Eqn. (A5). We can define an average replication rate for 
the single strand as 

A = (A6) 

yielding the main result of this appendix, 
dx^, 



dt 



= AmW{c^, cj)')x^, - f{t)x^. (A7) 



which looks remarkably similar to Eqn. (Al). For large 
populations, A{(j)') must rapidly equilibrate and remain 
steady, yielding a system of equations that are identical 
to that of a conservative quasispecies with a transformed 
set of rate constants. This is studied in more detail in 
the particular example of appendix B. 



APPENDIX B: SINGLE FITNESS PEAK 
LANDSCAPE 

In this appendix, we study the system of appendix 
A evolving on the commonly used single fitness peak 
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landscape described below. To solve this problem, we 
shall explicitly make a number of approximations that 
have been well studied and found to accurately describe 
the true dynamics of the system for reasonable genome 
lengths and population sizes. 

The single fitness peak landscape describes the situ- 
ation where a specific genome perfectly "fits" the en- 
vironment and hence replicates rapidly, while all other 
genomes are equally poor replicators. Here, we investi- 
gate the case where at least one strand needs to be per- 
fect in order to be viable, a reasonable model for the 
housekeeping genes responsible for cell survival (other 
landscapes will be studied in a future work^). Thus, 
■^i't>a,4'b) = cr 1 if either (f>a or (f>b are in the set 
{4>o, (j^'o} which represent the master sequence and its per- 
fect complement, and A((j)a,4'b) — 1 otherwise. 

For large populations and genome lengths, we can ig- 
nore mutations from unfit sequences to the master se- 
quence (an approximation that becomes exact as the 
genome length increases to infinity, but is accurate at 
realistic finite genome lengths) and assume that, at equi- 
librium, master genomes are paired with statistically dis- 
tributed complements (which is exact for this landscape 
in the large population limit, and rapidly converges for fi- 
nite populations). Thus, using Eqns. (A6) and (A7) and 
the fact that the symmetric equations conserve the equal- 
ity of concentrations of complementary sequences, we 
can write a diff'erential equation for the sum of the pop- 
ulation fractions of the single stranded master genome 
and its complement, xq, and the remaining population, 
xi = 1 - xo, 

dxo 



when Xo = in the quasispecies solution. This gives 



-J- =q oxq~ fit)xo 
at 



dxi n , 



(Bl) 



Xo,l 



-f{t)xi 



(B2) 

where N is the length of the genome, q represents the 
replicative fidelity, or 1 — e, where e is the per base point 
mutation probability, which is assumed to be sequence 
independent, xq^i represents the fraction of the popu- 
lation that are imperfect sequences bonded to a per- 
fect sequence (that is, the strands that are not mem- 
bers of xq, but are bonded to a member of xo), and 
f{t) = axo + axo^i + 1 — xq — xo,i- As complements 
are statistically distributed, we can define xq,i — Fdixg, 
where i^o,i represents the fraction of perfect sequences 
bonded to imperfect sequences, and is independent of 
Xq. We can solve these equations by searching for equi- 
librium solutions, Xq = xi = 0, yielding two solutions, 
the quasispecies solution. 



1 



= (l + Fo,0(a-l) 

and the quasispecies-free solution, — 0. The error 
catastrophe occurs when the two solutions meet, i.e.. 



llq\ 



(B4) 



which is identical to the conservative solution and is plot- 
ted and discussed in section IL 



APPENDIX C: DEGRADATION OF 
POST-METHYLATION DNA REPAIR 

In this appendix, we evaluate the effect of the pro- 
gressive failure of post-methylation DNA repair on the 
single fitness peak landscape discussed in the previous 
appendix. Setting the probability that an error will be 
repaired to be A, we get 



E(A(0, (l)a)x{^^^^} -f A{(f)a, (j))X{^^^^}), (CI) 



where W(0, (^a,A) represents the A-dependent probabil- 
ity that unzipped strand (f)a will produce new strand 0, 
while the new quantity W2 (0, 4>a , A) represents the prob- 
ability that, after replication and post-methylation lesion 
repair of unzipped strand 4>a, the erroneous repair of er- 
rors will change 0a to strand cj). A set of manipulations 
similar to those in appendix A and the definition in Eqn. 
(A6) can be used to yield the equations 



dxj^ 
~dt' 



J2 AmiWicI^, 0', A) + W2{<l>, 0', X))x^' 



ifit)+A{(^))x^ 



(C2) 



When appfied to a single fitness peak landscape, these 
equations can be written as 



dxo ^ ^.N 

dt 2' 



-{f{t)+a)xo 



(1 - (1 - ^)^)^}^^o 



Xl 



1- Xq 



(C3) 



(C4) 



where xq and xi and f(t) are as defined in appendix B. 
This can be solved for the error threshold, which occurs 
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when approaches the full semiconservative treatment and the 

solution of appendix B in the limits A — > 1 and A — > 0, 



1 



a = 



(1 - Ae/2)^ + (!-(!- A/2)e)^ - 1 " 
This expression is plotted and discussed in section II, and 



(C5) respectively. 
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